Expanding universe as a classical solution in the Lorentzian matrix model 
for nonperturbative superstring theory 
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Recently we have shown by Monte Carlo simulation that expanding (3+l)-dimensional universe 
appears dynamically from a Lorentzian matrix model for type TIB superstring theory in (9+1)- 
dimensions. The mechanism for the spontaneous breaking of rotational symmetry relies crucially 
on the noncommutative nature of the space. Here we study the classical equations of motion as 
a complementary approach. In particular, we find a unique class of SO(3) symmetric solutions, 
which exhibits the time-dependence compatible with the expanding universe. The space-space 
noncommutativity is exactly zero, whereas the space-time noncommutativity becomes significant 
only towards the end of the expansion. We interpret the Monte Carlo results and the classical 
solution as describing the behavior of the model at earlier time and at later time, respectively. 

PACS numbers: 11.25.-w; 11.25.Sq 



Introduction. — It is widely believed that the birth 
of our universe can be described by superstring theory, 
which is a natural candidate for a unified theory includ- 
ing quantum gravity. Indeed, a lot of insights into this is- 
sue have been obtained by string cosmology over the last 
decade These studies are based on perturbative for- 
mulations incorporating nonperturbative effects through 
D-branes. An obvious drawback in such an approach, 
however, is that one has to choose a particular string 
vacuum from numerous vacua that are theoretically al- 
lowed. On the other hand, there is also a possibility 
that one can actually determine the true string vacuum 
uniquely if one uses a nonperturbative formulation. 

Along this line of thought, we have studied a S0(9,l) 
symmetric Lorentzian matrix model, which is considered 
to be a nonperturbative definition of type IIB superstring 
theory in (9-1-1) dimensions Q. Surprisingly our Monte 
Carlo results provide clear evidence that 3 out of 9 direc- 
tions start to expand at some critical time. The observed 
spontaneous breaking of the S0(9) rotational symme- 
try down to SO (3) has been understood intuitively by a 
mechanism, which relies crucially on the noncommutative 
nature of the space. While this is certainly intriguing, it 
also poses a crucial question whether the space-time be- 
comes commutative at later time as we observe it now. 

In this Letter we study the classical equations of mo- 
tion of the model as a complementary approach. In par- 
ticular, we find a unique class of S0(3) symmetric so- 
lutions, which turns out to have the time dependence 
compatible with the expanding universe. For this solu- 
tion, the space-space noncommutativity is exactly zero, 
whereas the space-time noncommutativity becomes sig- 
nificant only towards the end of the expansion. 

Lorentzian matrix model. — The matrix model pro- 
posed as a nonperturbative formulation of type IIB su- 



perstring theory has the action S = Sh + Sf, where [3] 
Si = -^tr[^a{CT'^)^p[A^,^fi]) , (1) 



with Afj^ (/^ 



= 0, • • • , 9) and (a = 1, • • • , 16) be- 
ing N X N traceless Hermitian matrices. The Lorentz 
indices /i and v are contracted using the metric rj = 
diag(— 1, 1, • • • ,1). The 16 x 16 matrices are ten- 
dimensional gamma matrices after the Weyl projection, 
and the unitary matrix C is the charge conjugation ma- 
trix. The action has manifest S0(9,l) symmetry, where 
Af^ and transform as a vector and a Majorana-Weyl 
spinor, respectively. The space-time is represented dy- 
namically by the ten bosonic matrices ^. 

An important feature of the Lorentzian model is that 
the bosonic part of the action is proportional to 



tr {F^.F'^'') = -2tr {Fo^f + tr (F,,)' 



(2) 



where F^i, = — z[Ap,v4,y] are Hermitian matrices, and 
hence the two terms in ^ have opposite signs. A com- 
mon approach to study the nonperturbative dynamics of 
this model was to make the Wick rotation Aq — iAio 
and to study the SO(IO) symmetric Euclidean model, 
which is proved to have finite partition function 0] ■ 
See Ref. [8| and references therein for studies of the spon- 
taneous symmetry breaking (SSB) of the SO(IO) in the 
Euclidean model On the other hand, it is suggested 
that the Lorentzian signature of the metric plays an im- 
portant role in the dynamics of quantum gravity 1^, 113 1 . 

In Ref. 0] we studied, for the first time, the nonper- 
turbative dynamics of the Lorentzian model defined by 



Z = j dAd^ e'^ = j dAe'^'"ViM{A) 



(3) 
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where the PfafBan PfyVl [A] appears from integrating out 
the fermionic matrices V^q. We made the partition func- 
tion ([3]) finite by introducing infrared cutoffs in both the 
spatial and temporal directions instead of making the 
Wick rotation. It was shown by Monte Carlo simulation 
that one can remove these cutoffs in the large- TV limit, 
and that the theory thus obtained has no parameters 
other than one scale parameter. 

The classical solutions. — Taking account of the in- 
frared cutoffs introduced in the Lorentzian model, we 
search for stationary points of the bosonic action Sh for 
fixed jftr (Ao)"^ and ■^tr(Ai)^. Then the problem re- 
duces to solving the classical equations of motion 

- [Ao, [Ao,A,]] + [A,, [A,, A,]] - AA, = , 

-[A,,[Aj,Ao]]-~XAo - , (4) 

where A and A represent the Lagrange multipliers cor- 
responding to the constraints. We look for solutions, 
which are given by a unitary representation of a Lie al- 
gebra [A^jAj/] = iffxuxAx, which guarantees automati- 
cally that the Jacobi identity is satisfied. (See Ref. [14] 
for an analogous study in the Euclidean model.) Moti- 
vated by the Monte Carlo results mentioned above, we 
restrict ourselves to solutions with ^/ = (4 < / < 9) 
and with SO (3) symmetry corresponding to rotations in 
the i — 1,2,3 directions. From the complete list of real 
Lie algebras with four generators (0 < /i < 3), the one 
with SO (3) symmetry is given uniquely by [eo, e^] = —i 
for i = 1,2,3, whereas all the other commutators van- 
ish. (This corresponds to the algebra Al\ in Table I of 
Ref. [il] for a = = 1.) 

The unitary irreducible representations of the above 
algebra are classified into two categories. One consists 
of the trivial one-dimensional representations given by 
eo = a and et — 0, where a is a real parameter. The 
other consists of the infinite-dimensional representations 
given by the operators cq — —i^ and e.^ — Oi exp(a;) on 
the space of functions of x with integrability, where 
the three real parameters ai specify a representation. As 
the basis of the functional space, we use the eigenfunc- 
tions of the Hamiltonian of a one-dimensional harmonic 
oscillator, which are given as 

U{x) = cMx) e-^-' , c„ = (^i/4v^2"/2)-i . 

The representation matrices of cq and ei/oi, which we 
denote as P and K, respectively, have the following ele- 
ments. 

Pnm = / dxfn(x)*{-i)-^f,n{x) 



^nni j dxfrfii^X^ e friii,X) 

= CnCme^/'^ J dxe^"" Hn{x+^)Hm{x + ^) 

= e^^^2-^''-"'^/^V^. 

M _i 
x^[2'/!(Af-0!(ki-™|+0! 

where M = min(n,m). In the last equality, we have 
used the property ^_ff„(x) = 2nHn-\ of the Hermite 
polynomials. 

Using a direct sum of the non-trivial representations, 
we find a set of SO (3) symmetric solutions to (|4]), which 
is given by 

Ao = \/AP(8)lfe , (5) 
Ai = K (g) diag(xH, • • • , xfcj) . (6) 

The parameters Xai = (xa)i should be chosen such that 
the points (a = 1, • • • , fc) have spherically symmet- 
ric distribution in the 3-dimensional space. One of the 
Lagrange multiplier is fixed as A = 0. 

In the following analysis, the k x k matrices that ap- 
pear in ([5]) and ([6]) are omitted since they only give an 
irrelevant constant factor. Also we consider only one spa- 
tial direction z = 1 for simplicity since it turns out that 
the number of spatial directions does not play any role. 

The space-time structure. — In order to extract the 
space-time structure from the solution, we first need to 
diagonalize Aq. We do this numerically by truncating the 
functional space to the A^-dimensional space spanned by 
fn{x) with < n < — 1. Let us define the eigen- 
vectors \ta) corresponding to the eigenvalues ta of Aq 
(a = 1, • • • , A^) with the specific order ti < ■ ■ ■ < t^- 
The spatial matrix {ti\Ai\tj) in that basis is not diago- 
nal. However, it turns out that the off-diagonal elements 
decay exponentially in the direction orthogonal to the 
diagonal line. 

To see it explicitly, let us consider the N x N matrix 
QiJ = {ti\{Ai)'^\tj). Plotting yiQ7j7Qjv/vv/2l against 
I— J for a fixed value of (/+ J)/2, we find that it decreases 
exponentially with |/ — J|. The half width is largest 
for {I + J)/2 = N/2, and we denote it as n for later 
convenience. For A^ — 16, 32, 64, 128, we obtain n — 
11,15,23,33. 

The above observation motivates us to define n x n 
matrices A[°'''\t) = {t,y-f-a\Ai\ti,^b) with 1 < a,b < n and 
t = ^ J2a=i ^v+a iov V — Q, - ■ ■ ,{N — n). These matrices 
represent the space structure at fixed time t. Let us define 
the extent of space at the time t as R{t)'^ = ^ti Ai{t)'^. 
In Fig. [U we plot i?(i)/i?(0) for A^ = 16, 32, 64"l28. It is 
symmetric under the time reflection t ^ —t as one can 
prove analytically even at finite A^. For each A^, we have 
chosen the Lagrange multiplier A, which determines the 
scale of t, so that R{t) scales around t = 0. We have fixed 
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A = 1 for A'^ = 16 without loss of generality. Then we ob- 
tain A = 0.92, 0.72, 0.59 for iV = 32, 64, 128, respectively. 
As we increase N, the scaling region extends to larger 
\t\. The solid line is a fit to the Gaussian function. Thus 
we find that the time-evolution of the space is compati- 
ble with the expanding behavior observed in the Monte 
Carlo simulation Q. 
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FIG. 2: The dimensionless parameter representing the 
space-time noncommutativity is plotted against \t\ for A*' = 
16, 32, 64, 128. We have used the same set of values of n and 
A for each A'^ as in the previous figure. The data show nice 
scaling in A'^, and the solid fine represents x = 0.629 
obtained by fitting the A'' = 128 data within —6 < t < —2 to 
a power law. 



FIG. 1: The extent of space R{t) / R{0) is plotted as a function 
of t for four values of A'. The block size n is determined from 
the the decay rate of the off-diagonal elements of Ai in the 
basis which diagonalizes Aq . The value of A is chosen for each 
A' in such a way that the results scale in A'^. The solid line 
represents y = exp(— 0.034 f^), which is obtained by fitting 
the A' — 128 data to the Gaussian function. 

Let us next turn our attention to the space-time non- 
commutativity. We define the dimensionless parameter 



(7) 



as an estimate on the space-time noncommutativity [16| . 
In Fig. m we plot for N = 16,32,64,128. We find 
that it is of 0(1) at t = and decreases as ~ at large 
\t\, where we obtain p = 1.708(3) by fitting the = 128 
data within —6 < t < —2. Therefore, the space-time 
noncommutativity is significant only around t '-^ 0, and 
it becomes smaller as we go back in time. 

In both Figs. [1] and O we consider it very important 
to study the larger |t| region by increasing N. In fact 
some preliminary results suggest certain deviation from 
the Gaussian behavior of R{t) and the power-law behav- 
ior of x(i)- We hope to report on it in future publications. 

Summary and discussions. — We have studied the 
classical equations of motion in the Lorentzian matrix 
model for type IIB superstring theory. Restricting our- 
selves to the class of solutions that are written in terms 
of Lie algebras with four generators, we find that the so- 
lution with SO (3) symmetry is essentially unique. The 
space-time structure extracted from the solution exhibits 



the time dependence, which is compatible with the be- 
havior observed in our previous Monte Carlo results. The 
space-time noncommutativity becomes significant only 
towards the end of the expansion, whereas the space- 
space noncommutativity is zero. These results suggest 
the appearance of an expanding (3-|-l)-dimensional (al- 
most commutative) space-time from the Lorentzian ma- 
trix model. 

We speculate that the noncommutativity of space 
which plays a crucial role in making three directions ex- 
pand at earlier time, somehow disappears at some point 
due to some dynamical reason. For instance, let us use 
the model obtained after integrating out the scale factor 
Q . In that model we have a constraint that requires the 
quantity ([2]) to vanish. If the expansion with large space- 
space noncommutativity in the early time continues for 
too long a period, the second term of ^ will be too large 
to satisfy the constraint ^= 0. Such an effect may lead 
to an end of the noncommutative expansion. One might 
speculate that this corresponds to the end of "inflation" . 

Our classical solution is symmetric under time reflec- 
tion, and the size of the space becomes maximum at 
t = 0, after which it has a contracting behavior. At i = 0, 
the dimensionless space-time noncommutativity becomes 
maximum, too, and it is of the order one. Hence the 
physics there will be quite exotic. This may be taken as a 
prediction on the fate of our universe from the Lorentzian 
matrix model given that our classical solution is valid 
around t = 0. 

Obviously one can generalize our solution to SO(c?) 
symmetric ones with 1 < c? < 9. The time evolution 
of the size of the space and that of the space-time non- 
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commutativity are essentially the same as in the S0(3) 
case. Let us recall here that the space-space noncommu- 
tativity seems to play a crucial role in the SSB of S0(9) 
rotational symmetry as the mechanism proposed in our 
previous work Q suggests. It is therefore not surprising 
that the dimensionality of space is not fixed by classical 
solutions without space-space noncommutativity. 

Some comments on related works are in order. In 
Ref. the Einstein equation was derived from the clas- 
sical equations of motion of the matrix model. In this 
derivation, however, the matrices were interpreted as 
the covariant derivative on a curved space. It would be 
interesting to clarify the relationship to our work. 

Ref. ^9(1 reports on interesting solutions to the classical 
equations of motions (|4]) with A = A = 0. They represent 
a flat Minkowski space with extra dimensions described 
by fuzzy spheres. An interesting feature of these solutions 
is that there exists noncommutativity between the four 
extended directions and the extra dimensions. This is 
crucial for realizing a nontrivial structure in the extra 
dimensions even without the Myers-like term. Here we 
emphasize that the nonzero A, which is introduced in our 
work, is crucial for the expanding behavior. Let us recall 
that A is the Lagrange multiplier corresponding to the 
infrared cutoff, which turns out to be needed according 
to our previous Monte Carlo studies. 

In Ref. [i^l the Matrix theory 21 [ has been applied 
to cosmology. A classical solution with three expanding 
(commutative) directions and six oscillating (noncommu- 
tative) directions was discussed. (The number of expand- 
ing directions does not have to be three.) In order to have 
such a solution, the authors introduced a SO (9) symmet- 
ric tachyonic mass term, which was interpreted as the 
cosmological term. The relationship to our solution is 
not clear, though, since the time is treated in a different 
way. The idea to use the matrices to avoid the big-bang 
singularity is also pursued in Refs. (22| . 

It is tempting to imagine that the rapid growth of R(t) 
observed in the present solution has something to do with 
the accelerating expansion confirmed by recent cosmo- 
logical observations. The power-law expansion at earlier 
time may be understood by considering the quantum cor- 
rections around the classical solution. It is also expected 
that the gauge interactions and the matter content in the 
(3-|-l)-dimensional space are determined by the structure 
in the extra dimensions 231 l2J| analogously to the case of 
intersecting D-brane models. We hope the present model 
provides a new perspective on particle physics beyond the 
standard model as well as on cosmological models for in- 
flation, modified gravity etc.. 
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